Three-body direct nuclear reactions: Nonlocal optical potential 
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The calculations of three-body direct nuclear reactions with nonlocal optical potentials are per- 
formed for the first time using the framework of Faddeev-type scattering equations. Important 
nonlocality effect is found for transfer reactions like d + ^^O — > p-|- ^'^O often improving the descrip- 
tion of the experimental data. 

PACS numbers: 24.10.-i, 21.45.-v, 25.55.Ci, 24.70. -|-s 
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Direct nuclear reactions dominated by three-body de- 
grees of freedom provide an important test for mod- 
els of nuclear dynamics. Extensively studied examples 
are deuteron [d) scattering from a stable nucleus {A) 
and proton {p) scattering from a weakly bound nucleus 
{An) consisting of a core A and a neutron [n). The 
nucleon-nucleus {NA) interactions employed in three- 
body calculations are usually modeled by optical po- 
tentials (OP) with local central and spin-orbit parts. 
This local approximation yields a tremendous simpli- 
fication in the practical realization of the continuum 
discretized coupled channels (CDCC) method the 
distorted- wave Born approximation (DWBA), or vari- 
ous adiabatic approaches [1| that are widely used for the 
description of three-body nuclear reactions; to the best 
of our knowledge, the above mentioned methods using 
nonlocal optical potentials (NLOP) have never been at- 
tempted. However, NLOP can be included quite easily in 
the framework of exact three-body Faddeev/Alt, Grass- 
berger, and Sandhas (AGS) scattering theory if one 
uses the momentum-space representation. Its application 
to three-body nuclear reactions [5, 6] has become possible 
recently due to a novel implementation of the screening 
and renormalization method 0, H, @ for including the 
long-range Coulomb force between charged particles. 

The aim of the present work is to calculate the ob- 
servables oi d + A and p + {An) reactions in a three- 
body model (n,p, A) using NLOP and its (nearly) equiv- 
alent local optical potential (LOP) in the framework of 
momentum-space AGS equations and thereby study the 
effect of the OP nonlocality. 

The AGS equations [1] are Faddeev-like connected- 
kernel equations that provide an exact description of the 
quantum three-body scattering problem. In contrast to 
the Faddeev equations Q formulated for the components 
of the wave function, the AGS equations, 
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(1) 

are a system of coupled integral equations for the tran- 
sition operators Upa{Z) whose on-shell matrix elements 
are scattering amplitudes and therefore lead directly to 
the observables. In Eq. ([T]) S(3a = 1 — 6pa, Go{Z) = 
{Z — Ho)~^ is the free resolvent, and 
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is the two-particle transition matrix with Z = E + iO, 
E being the available three-particle energy in the cen- 
ter of mass (cm.) system, Hq the free Hamiltonian, and 
v-y the potential for the pair 7 in odd-man-out notation. 
We work in momentum space where local and nonlocal 
potentials are treated in the same manner. The AGS 
equations are solved numerically after partial-wave de- 
composition and discretization of momentum variables; 
more details on ernployed numerical techniques can be 
found in Refs. [13, [ill- 

The AGS equations were formulated originally for a 
short-range potentials v-,. Nevertheless, the method of 
screening and renormalization 0, H, [§] allows to include 
the long-range Coulomb force between charged particles 
using the AGS framework. The Coulomb-distorted short- 
range part of the transition amplitude is obtained by 
solving the AGS equations with nuclear plus screened 
Coulomb potentials; the convergence of the results with 
the screening radius has to be established. The method 
has been successfully applied to proton-deuteron [H, US] 
elastic scattering and breakup and to three-body nuclear 
reactions involving deuterons or one-neutron halo nuclei 

m- 

For the np interaction we take the realistic CD Bonn 
potential |14| , in contrast to a simple Gaussian np poten- 
tial used in CDCC, DWBA, and adiabatic calculations. 
For the hadronic part of the NA interaction we take the 
NLOP of Giannini and Ricco 15] which in configuration 
space has the form 

«^(r',r) - H{x)[V{y) + tW{y)] + cr ■ \.R,{x)V,{y) (3) 

with X — |r' — r|, y = |r' -|- r|/2, 
H{x) = (7r/32)-3/2cxp(-a;V/32), H,{x) 

(^/32)-3/2exp(-a;V/3s'); "^{v). and V,{y) 

are the real volume, imaginary surface, and real spin- 
orbit parts, respectively, that are parametrized in the 
standard way using Woods-Saxon functions. The ap- 
proximately equivalent LOP is taken over from Ref. [15| 
as well. Both NLOP and LOP are slightly modified: 
we adjust the parameter Wjv, determining the strength 
of the imaginary part, to improve the description of 
the experimental N A scattering data in the considered 
nucleon lab energy range from 10 to 40 MeV as well as 
the agreement between the N A predictions of NLOP 
and LOP. The adjusted values for Wj^ are given in 
Table HI the other parameters are taken from Ref. [lH| ■ 



2 



TABLE I: The parameter Wn of NLOP and LOP (in units 
of MeV) adjusted to the NA data. 





Wiv(NLOP) 


Wn (LOP) 




13.0 


10.0 




14.0 


11.0 




15.0 


12.0 



In contrast to the NLOP, the LOP is energy-dependent 
owing to the equivalence transformation [l^. Though 
energy-dependent potentials were used recently [T^ in 
Faddeev-type calculations, we refrain to do so in the 
present work. We follow the standard procedure of 
fixing the energy-dependent parameters of the two-body 
OP in three-body calculations and use two types of 
Hamiltonians: (a) In d -I- ^ elastic scattering both nA 
and pA OP parameters are taken at half deuteron lab 
energy, (b) In p + {An) elastic scattering and transfer 
to d -I- ^ the parameters of the pA OP are taken at the 
proton lab energy whereas the nA potential has to be 
real in order to support an (An) bound state. Since the 
observables in p + (An) reactions are rather insensitive 
to the nA potential, provided it reproduces the spectrum 
of bound states, we take over the local real nA potential 
from Ref. [l3| that supports a number of bound states 
corresponding to the ground and excited single-particle 
states of the (An) nucleus, while all Pauli-forbidden 
states are removed; the potential parameters and the 
resulting binding energies are given in ReL [l^ for 
^^C and ^"^O nuclei. The standard Hamiltonian (a) 
with the nA potential being complex does not support 
{An) bound states and therefore does not allow for 
calculations oi d + A ^ p + (An) reactions while most 
of the available transfer data come from this type of 
reactions. However, since the d + A ^ p + (An) and 
p + (An) d + A reactions are related by time reversal 
provided the energy in the cm. system is the same, we 
calculate the latter one using the standard Hamiltonian 
(b) where the nucleus (An) can be in its ground or 
excited state and apply the time reversal to obtain the 
observables for the former one. This is equivalent to 
using the Hamiltonian (b) in the d + A scattering, a 
nonstandard choice. 

The interaction between np, nA, and pA pairs is in- 
cluded in partial waves with pair orbital angular momen- 
tum L < 3, 10, and 20, respectively, and the total angular 
momentum is J < 40; depending on the reaction some of 
these quantum numbers cutoffs can be safely chosen sig- 
nificantly lower, leading, nevertheless, to well converged 
results. The pA channel is more demanding than the nA 
channel due to the screened Coulomb force, where the 
screening radius i? « 8 to 10 fm for the short-range part 
of the scattering amplitude is sufficient for convergence. 

In Fig. [T] we use p + ^^O elastic scattering at proton 
lab energy Ep = 35.2 MeV as an example to illustrate 
the achieved quality in fitting the two-body data and the 
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FIG. 1: (Color online) Differential cross section divided by 
Rutherford cross section for p-|-^^0 elastic scattering at Ep = 
35.2 MeV. Predictions of NLOP (solid curve), LOP (dashed- 
dotted curve), and Watson OP (dotted curve) are compared 
with the experimental data from Ref. [l7l |. 
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FIG. 2: (Color online) Differential cross section divided by 
Rutherford cross section for p + ^^O elastic scattering at 
Ep = 35.2 MeV. Predictions of NLOP (solid curve) and LOP 
(dashed-dotted curve) are compared with the experimental 
data from Ref. |17| . 

approximate NLOP-LOP equivalence. We show also the 
predictions of the OP by Watson et al. [18] to demon- 
strate that the adjusted NLOP and LOP describe the 
NA data at least as good as the traditionally used po- 
tentials. An agreement of similar quality is found in all 
considered cases. 

In Figs. [2] and [3] we show the results of the three-body 
calculation for elastic proton scattering from ^^O and ^^C 
nuclei around Ep = 35 MeV. The mutual agreement be- 
tween NLOP and LOP predictions and with the data is 
as good as in the two-body case shown in Fig. [TJ This is 
not surprising since the differential cross section and the 
proton analyzing power in p + {An) elastic scattering are 
known to correlate strongly with the corresponding ob- 
servables \np + A elastic scattering. Therefore, the effect 
on the OP nonlocality is tiny; even the small differences 
between NLOP and LOP predictions seen in p+^^0 cross 
section at large angles are reproduced well in p-f ^^O elas- 
tic scattering. 

The Hamiltonian (b) used to calculate p -\- ^^C elas- 
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FIG. 3: (Color online) Differential cross section divided by 
Rutherford cross section and proton analyzing power for for 
p+^^C elastic scattering at Ep = 35 MeV. Curves as in Fig. [2] 
and the experimental data are from Ref. [l9l |. 
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FIG. 4: (Color online) Differential cross section for p + ^^C — > 
d+^'^C transfer at Ep = 35 MeV. Curves as in Fig. [5] and the 
experimental data are from Ref. [20| ]. 



tic scattering in Fig. [3] allows for the transfer to d + ^^C 
as well; the AGS equations for both reactions are solved 
simultaneously. In contrast to elastic scattering, the ef- 
fect of the OP nonlocality becomes significant in the 
p+^^C d+^^C transfer cross section at angles Qc.m. > 
30 deg as demonstrated in Fig. [4l although it does not 
improve the description of the experimental data. We 
therefore expect the OP nonlocality to be important also 
in the inverse reactions, i.e., d + A p + {An), where 
more data exist with the final nucleus (An) being in its 
ground or excited state; the observables are calculated 
using the Hamiltonian (b) as described above. We start 
in Fig. [5] with the d+^^C — > p+^^C reaction at deuteron 
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FIG. 5: (Color online) Differential cross section for d+^^C —* 
p+^'^C transfer to ground and excited '^^C states at Ed = 30 
MeV. Curves as in Fig. [2] and the data are from Ref. [2ll ]. 



lab energy Ed = 30 MeV which in the p + ^^C ^ d + 
case corresponds to Ep — 30.6, 27.3, and 26.5 MeV for 
states 1/2", 1/2+, and 5/2+, respectively. For the 
transfer to the ^^C ground state 1/2^ the shape of the 
experimental data and theoretical predictions are similar 
to the ones of p+^^C d+^^C reaction at Ep = 35 MeV 
in Fig. m as expected from the detailed balance, given the 
small difference in energy. Thus, except for forward an- 
gles, the transfer reactions involving the ^^C ground state 
1/2" are described rather unsuccessfully much like it was 
found in Ref. In the case of the transfer to ^'^C ex- 
cited states 1/2+ and 5/2+ the OP nonlocality is again 
important, and the predictions of the NLOP account for 
the data quite successfully, in contrast to the ones of the 
LOP and Ref. [l^. The differential cross section for the 
transfer to "'^^C 1/2+ state is increased by the NLOP at 
forward angles and decreased at Qc.m. > 20 deg while for 
the 5/2+ state it is significantly decreased at Qc.m. > 35 
deg such that the data is slightly overestimated by al- 
most a constant factor that may be associated with the 
spectroscopic factor. 

The results for the d+^^0 p+^'^O transfer cross sec- 
tions at Ed = 25 A and 36.0 MeV are presented in Figs.jH] 
and [71 The corresponding proton lab energy in the in- 
verse reactions p + ^'^O ^ d + ^"^O is Ep = 25.9 and 35.9 
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FIG. 6: (Color online) Differential cross section for d+^^O —^ 
p + ^'^O transfer to the ground state 5/2+ of at Ed = 25.4 
and 36.0 MeV. Curves as in Fig. [2l The experimental data 
are from Ref. [H. 
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FIG. 7: (Color online) Differential cross section for d+^^O — > 
p + transfer to the excited state 1/2+ of ^''O at Ed = 25.4 
and 36.0 MeV. Curves as in Fig. [2l The experimental data 
are from Ref. [H. 



MeV for ground state 5/2+ and Ep = 25.0 and 35.0 
MeV for ^^O excited state 1/2+ , respectively. The effect 
of the OP nonlocality is again sizable and qualitatively 
very similar to the one observed in d + ^^C ^ p + ^^C 
reactions involving ^'^C excited states 5/2+ and 1/2+ . 
The differential cross section for the transfer to the ^^O 
ground state 5/2+ is decreased at very forward angles 
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FIG. 8: Differential cross section divided by Rutherford cross 
section for elastic scattering at Ed = 25.4 MeV. Re- 

sults obtained with the Hamiltonian of type (a) and (b) are 
shown on the left and right side, respectivel y. Curves as in 
Fig. [2] The experimental data are from Ref. |22l ]. 
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FIG. 9: Differential cross section divided by Rutherford cross 
section and deuteron vector analyzing power for deuteron 
elastic scattering from ^^O and ''"Ca nuclei at Ed = 56 MeV. 
Curves as in Fig. [5] The experimental data are from Ref. [23|. 



and at Qc.m. > 30 deg. Although in particular narrow 
angular regimes the predictions of the LOP are closer to 
the data, the NLOP provides a better overall description 
of the data, especially of its shape. In the case of the 
transfer to the ^^O excited state 1/2+ the NLOP repro- 
duces the data almost perfectly. 

Next we consider deuteron-nucleus elastic scattering 
that is usually described using the Hamiltonian of type 
(a). We show in Fig. [8] the results for d + ^^O elastic 
cross section at Ed = 25.4 MeV. The effect of the OP 
nonlocality is visible at larger angles Oc.m. > 30 deg and 
moves the predictions towards the data up to ©cm. — 60 
deg; beyond 60 deg the data is strongly underpredicted 
by both NLOP and LOP results. This indicates that the 
Hamiltonian (a) is too absorptive in the considered case 
where negative energies in the two-body NA subsystem 
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have a large weight in the three-body scattering equa- 
tions as discussed in Ref. [lB|; it was found there that 
using energy-dependent potential that becomes real and 
therefore less absorptive at negative NA energies signifi- 
cantly increases the cross section at 9c. m. > 60 deg. For 
curiosity in Fig. [8] we present also results of the Hamilto- 
nian (b) which is less absorptive since the nA potential 
is real. Surprisingly, with this choice both the NLOP 
and the LOP roughly account for the data also at large 
angles. Thus, NLOP and LOP with the Hamiltonian of 
type (b) describe elastic d + ^^O and p + ^^O data with 
similar quality while for the transfer reactions the NLOP 
is clearly superior. 

Finally in Fig. [5] we present observables for d + ^^O 
and d + ''"Ca elastic scattering at Ed = 56 MeV calcu- 
lated with the Hamiltonian of type (a). Visible differ- 
ences between NLOP and LOP show up at large angles 
©cm. > 60 deg. The predictions of NLOP are closer to 
the data but show too sharp oscillations; a similar feature 
is present already in the p + A and p + (An) observables 
in Figs.[T] — [3]but is less pronounced there. Both NLOP 
and LOP fail in reproducing small-angle deuteron vector 
analyzing power Ay much like other potentials as found 
in Ref. ^24| . The effect of the OP nonlocality is of compa- 
rable size also for deuteron tensor analyzing powers, but 
it is unable to resolve the discrepancy [13] in the large- 
angle Axz- The results for deuteron-^^C elastic scattering 
exhibit similar features as those for deuteron- ^^O and are 
not shown separately. 

To be sure that the observed OP nonlocality effect 
is not a consequence of only approximate NLOP-LOP 
equivalence, we performed the following test calculations. 
We varied the parameter Wn by ±1 MeV thereby induc- 
ing changes in the NA observables of a size comparable 



to the NLOP-LOP difference. As a consequence, sim- 
ilar changes occur in p + (An) elastic observables, but 
for d + A elastic scattering and especially for all transfer 
reactions the induced changes are considerably smaller 
than the observed NLOP-LOP difference. Thus, the im- 
perfection in the NLOP-LOP equivalence does not affect 
our conclusions on the importance of the OP nonlocality. 
A possible reason for an overall better description of the 
three-body observables by NLOP may be that NLOP fits 
the NA data over a broader energy range compared to 
LOP which, although being energy-dependent, has to be 
chosen at a fixed energy. 

In summary, we performed the calculations of three- 
body direct nuclear reactions with the NLOP for the first 
time. Exact scattering equations in the AGS form were 
solved and the Coulomb interaction was included using 
the method of screening and renormalization. The OP 
nonlocality effect is found to be very small for p + (An) 
elastic scattering, moderate for d + A elastic scattering 
at larger angles, and especially important for transfer 
reactions p + (An) d + A and d + A p + (An). 
In the latter case the NLOP is clearly more successful 
in accounting for the data in transfer reactions involving 
and ^^O nuclei in the states 1/2+ and 5/2+. We hope 
that the present work, demonstrating the feasibility of 
the calculations with NLOP and the importance of the 
nonlocality, will stimulate the development of new and 
more precise nonlocal nuclear interaction models. 
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